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Problem Description: 
To investigate the response of an automobile suspension system for selected disturbances. The 
system response to these disturbances for various values of the system design parameters are obtained. 
The most suitable values of the system parameters are determined by selecting the desirable response 
from the computer solution. 
a Where : 
M, = One quarter of mass of automobile 
M, = Mass of the wheel and axle 
K, = Spring constant of main auto spring 
K, = Spring constant of tire (assumed linear) 
D, = Shock absorber damping constant 
x, = Displacement of auto body 
x, = Displacement of wheel 
x; = Roadway profile displacement 
Physical Constants: 
M, = 25 slugs 
M, = 2 slugs 
K, = 1000 lbs/ft. 
K, = 4500 lbs/ft. 
D, = Variable 
ne Initial conditions and forcing functions: 
dx, dx. 
Figure 1: Simplified representation of My, = hy = = = 0 att =0 
a single wheel of an automobile sus- dt dt 
pension system. 
%, = x(t) 
System Equations: 
The differential equations of motion of the system are derived by equating the forces acting upon 
the mass involved in the system. 
They are: dx, Dy ( dx, dx, _ Ay a 
dt? my, dt dt Mm, 
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Scaled Equations : 


Assume the following maximum values of the scaled variables to assure that the computer problem 
voltages are as high as possible but less than + 10 volts. 





d2x% d2x 

© ee 0.8 : 

[ - ] [ aa" 
d 2 

[5 —*] [0.8 22] 
dt dt 

[ 50 x, | [ 50 x, | 


The scaled equations then become: 


[s%l=- (=) tsa)+ (7) [0.8%] - (4.0) [50%] + (4.0) [50x] 
25 4 
[0.8%,] = — — | [0.8 x,] + (0.08 D,) [5x,] — (44) [50 ~%,] + (8) [50 x,] + (1800) [x,] 
Where ze d2x a dx 
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Combination potentiometer settings and amplifier gains. 
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Voltage outputs from computing components. 


Computer Diagram: 





Fach input to every integrator is multiplied by a factor —to change the speed with which the 


system behavior is obtained on the computer. In this example B = 10 appears appropriate and the com- 
puter circuit behavior will be ten times slower than that of the suspension system. 
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Figure Ic. 


X Xx 
Displacement of auto body ify and wheel oe for a deflection in road profile (X, is a square pulse of 35 
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CALCULATION OF HEAT TRANSFER BY NATURAL CONVECTION: 


Problem Description: 

To solve a natural convection heat transfer problem, the solution will consist of the temperature and 
velocity distribution for various Prandtl numbers. The quantity of heat which is transferred will also be cal- 
culated. 


Equations : 


The Navier-Stokes equation and the Fourier equation can be written for natural convection heat 
ned transfer from a vertical flat plate in the following manner’. 


ay Oty oy | aM ay 


= ¥V 


dy Ox doy Ox dy? dy3 





+ gBO,¢ 


dy 6 _ ay a6 _ as 


dy Ox Ox dy dy? 


with boundary conditions 


; 
Woy Zeb B= 1 Me = B 








Oy Ox 
e.g 2s 5 oe - 0 At y = @ 
_ Oy Ox 
Using substitutions suggested by Pohlhausen! of : 
¢ = & é 
x A 
where 
(Sr 
Av2 
and le) = (x,y) 
W (x,y) 
2.) 
, Ave x% 


The original equations reduce to the following system of ordinary differential equations : 


4+3 48-2 272+ = 0 


7+ 3P, £7 = 0 
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With boundary conditions 


Ze = 0 Zed At « = 0 


II 
p— 


Z 0 At € = 


ll 

© 

tN: 

Il 

© 
3 

It 


where the dots denote differentiation with respect to e. 


In order to solve these equations a substitution is made for the independent variable. 


e = t/a a-z= 10 


dad a y an ot 
de * dt ™° de® ~ * dt2 


Computer Diagram: 





Results of Analog Computer Solution: 


The solution of the problem is one of trial and error in that all the initial conditions (the values of 
the functions ate = 0) must be known. Values for 4 and 7 ate = O are assumed and if the boundary con- 
ditions are satisfied, the assumed values are correct. 


Figures 1 and 2 show the temperature and velocity distributions for various values of the Prandtl 
number as computed. 


To calculate the quantity of heat which is transferred from the plate only the value of 7 ate = Ois 
required. The desired values of the design parameters of interest can then be calculated using well known 
steady state heat transfer relationships ! 
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: Heat transfer temperature profile. 
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Figure 2: Heat transfer velocity profile. 


Nomenclature: 


A T T 
g 8 Os ) p_ + +00 

( 4v2 2 
Cp = Heat capacity of fluid u Velocity in vertical direction 
8 B aL Vv 


Gp = Se (Grassoff number) 
Vv 





If 


ll 


Velocity in horizontal direction 


Vertical and horizontal coordinates of a 
point in fluid 

_k 

Cpp 


” 
“<j 
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g = Acceleration due to gravity 
h = Over-all heat transfer coefficient. o = 
k = Thermal conductivity of fluid 
L = Height of plate B = Poo ~ Pa 
Nu = ut (Nusselt number ) Palla = Te? 
wu = Dynamic viscosity 
Cpu = Xj ti . ‘+ 
P,, ~ ( Prandtl number ) Vv Inematic viscosity 
k = Density of fluid 
T = Temperature of fluid at point T — T,, 
Too = ‘Temperature of bulk fluid e = T T 
P ~~ “oo 
Tp = Temperature of plate 
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CALCULATION OF RADIAL VELOCITY OF ROTARY SPRAY DRIER 





Problem Description: 

To compute a curve of radial velocity of a particle vs radial distance from the center of a disc used 
in a spray drier. The results of the computation are then used to select the design parameters of a drier for 
a particular duty. 


Equations: 
lea The differential equation which describes the velocity of the particle is I. 
v *Y 4 av? _ Br = 0 
dr 
where 
V = Radial Velocity 
r = Radial distance from center of disc 
A, B = Constants which in the case of laminar flow depend upon the angular velocity of disc, volu- 
metric feed rate, vane height, and the density and viscosity of the liquid. For turbulent flow 
the same equation hold, with different value for constant A. 
Physical Constants: 
Mean values of 
A = 1.66 x 107+ 
B= 36 x 107 
r = Oto 5 inches 
V = 0 to 6000 ft/min 
NS 


Computer Equations: 


In order to solve this equation a substitution must be made for the independent variable r: 


let r = t/a 
then d d 
dr ~ dt 


The equation becomes with physical constants : 


dy . S610’. 166 x 1"%Vv" 
dt 5V 5 


where 
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The computer will compute 10~°V, thus the equation becomes 


lo° dv _10°Br — 10%Av? 
dt 5V 5 
Calculation of radius: 
l 
_ dr = —dt 
b= Fé a 2r = 2 fdt 
dt = adr a 
l 
r = — fdt 
a 


COMPUTER DIAGRAM 








. 


w 
a 


VELOCITY IN THOUSANDS FT/SEC 


VELOCITY IN THOUSANDS FT/SEC 


Nn 





seiaastteastaastositicaaiae 


3 * 5 . RADIUS JN INCHES 
RADIUS IN INCHES ra 


Figure 1: Plot of Radial Velocity as a Function of Drier Radius. 
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ANALYSIS OF THE FLOW PATH OF AN OIL GLOBULE 


Problem Description: 


It is required to separate an oil-water mixture by using a standard oil separator. The over-flowing 
level of oil is 5 feet above the lower end of the oil retention wall. Calculate the minimum length of separa- 


tion section of the separator. 
Equations : 


OIL—WATER To analyze the flow path of an oil globule, 
MIXTURE IN its velocity is resolved into vertical compo- 


( OIL RENTENTION SECTION nent v and horizontal component h. The hori- 






—> OVERFLOWING OIL zontal velocity of an oil globule is assumed 
constant (throughout the separation section) 
and equal to the total volume flow rate divided 
by the cross-sectional area of the separator. 
The vertical velocity is variable and depen- 


s—>WATER OUT dent upon a force balance which is derived as 









SS 








Fj 
' | v ie PATH 


Seri h 
a VOIL ~GLOBULE 






MMOOOOH}y 
RV] g$agy 
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ollows: 
SCHEMATIC OF STANDARD 
OIL—WATER SEPARATOR Inertial force of the oil globule 
D3 dv 
yoa- re 
se 2 
Buoyancy of the globule 
a D8 
Fb = : (Py — Po) 8 
Drag force against the upflowing globule 
aD? ‘ 
v 
Fa = d ( 4 de 
2 
Where 
D = Average diameter of oil globule 
v = Vertical velocity of oil globule 
po = Density of oil 
Pw = Density of water 
g = Acceleration due to gravity 
f = Coefficient of viscous friction 





Long Branch, New Jersey 





ELECTRONIC ASSOCIATES INC. 





een? BULLETIN #AN 924-1 


U.S. AF94i1 


The force balance becomes: 











F Fy = Fa 
a D3 dv a D3 f aD? py v? 
po — = ‘fo =~ pl sg - 
6 dt 6 8 
dv Pw — Po 3f Pw 
dt 7 Po 6 AD po 


The vertical velocity v may by integrated with respect to time to give the vertical traverse distance 


of an oil globule within a given time. 


t 
s=- ff @ dt 


0 

The computer diagram is mechanized to solve for v and s. Since the total vertical traverse distance 
is known and the function of vertical velocity will be obtained from the computer, the above equation is 
employed to determine the time required for an oil globule to rise to the over-flowing surface. The mini- 
mum longitudinal length of the separation section equals the horizontal velocity h times the minimum re- 
tention time. 

The variables are scaled as indicated below so that the computer problem voltages are as high as 
possible but less than + 10 volts. 


[100 iso = 
dt dt 
[ 100v ] [ 104s ] 


Scaled Equations: 


Pw — po  3f py LLOOV)? 


100 = ly Sic hc eke ence amie 
esi a _—" 40D po [101 


Computer Diagram: 
+ [104s ] 









10? 8 (Pw —Po) 
[ 100v]2 


. [10] 


#19 





3f Pay 
4D p Bx 10? 


200 


108 


= CO {> +2001 


] 
Each input to every integrator is multiplied by a time scale factor. In this example 8 = 100 and 


the computer solution will be 100 times slower than the behavior of the physical system. 


Computer Results : 
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Plot of vertical velocity of oil globule with time. 
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INVESTIGATION OF A SIMPLE CHEMICAL REACTION 


Problem Description: 


The reaction of chemical component A to form component B, which in turn reacts to form component 


C, can be represented as follows: 


A ry s B T, > 








where r, and r, are specific reaction rates. The reaction is assumed to take place at a constant tempera- 
ture. Find the changes in concentration of A, B and C as a reaction time, (for the case of r, equal to 


unity and r,/r,; less than 1.0.) 


Equations: 


The reaction rate for such reactions may be expressed by the following differential equations, 


assuming no reverse reactions: 





dn 
~ — = Thy rate of disappearance of A 
dn 
= yng — 7, np net rate of appearance of B 
dn 
— = Ty Bp rate of formation of C 
t 


Where n 4, np and n¢ are mole fractions. The mole fraction of a component is defined as the number of 


moles of that component divided’ by the sum of the number of moles of all components. 
Therefore O < nyg, np, ng < 1.0 
ng = 10 att = 0 
Scaled Equations: 


The variables are scaled as shown below to assure that the voltages representing their maximum 


values are as large as possible within the limits of + 10 volts. 


dn 4 dnp dnc 
[10 =] [10 —— | [10 —— 
[10n 4 | [10np | [10n¢| 


The scaled equations then become 


dn 4 
mz L10—— | = Ty [107 4 | 
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dnp 
[10 1] = Fy [lOey] —ry [10np] 


dnc 
[10 —*] = Ty, [10np | 


Computer Diagram: 


{tO 





+(1On,) 





+( lOng) 









Plot of Results Obtained: 
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Plot of concentrations as a function of reaction time for different specific reaction rates. 
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DYNAMIC RESPONSE OF A SEPARATELY EXCITED GENERATOR 


Problem Description: 


To investigate the time response of a particulardesign of a separately excited generator for various 


changes in the exciter voltage with a constant resistive load. 


Re Rf | Ra 


Te 





Schematic Diagram of Separately Excited Generator 


From the schematic the following equations can be written: 


=I,k. + Le — 
dt 


E.- 


L 


dl 
E. = Ip Re + Le as 


Eg = ly (Ra + Ry) 


Ry, 
E S same ie 
Ra + Ry, 


The armature voltage of the exciter EF, is proportional to the field current /,, so that 
fae = Beals 


Similarly the emf of the generator E, is proportional to the field current /f. 


Eg = Kg If 
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The equations to be programmed on the computer will then be: 


dE K R 

ne eee oe 

dt Le Le 

itp Ke, RY 

ag «(e Lf & Lf 

R 
6 means ff 
Ra + Rx 

Maximum values: E; = 200 volts 
E. = 200 volts 
Eg = 500 volts 
Ee = 500 volts 


Scaled Equations: 


The variables are scaled as shown below to assure that the voltages representing their maximum 
values do not exceed + 10 volts: 


1 dE 1 dE 
— —] [— —*] 
20 dt 50 dt 
E E 
[—] [] 
20 50 
E; Ey 
[— ] [ —] 
20 50 
and the scaled equations are of the form: 
1 dE, 
[— —*] = FS) (4 — GS (41 
20 dt 
1 dE 
ae EY te ey A mt 
50 dt 5 Lp Lf 
Eo R E 
i) = (dE 
50 Rag + Ry, 50 
Computer Diagram 
0 40-—+10 
+10 —10 











In+ [4] Be 
ie eed Eg Eo 
20 20 { [—] _ LS =. 
ua » <} 50 
wr HOO > >O->= 
Fi ot . “y RL 
200 5 Lf 





Results of Computer Solution: 


The response of the generator for changes in the exciter voltage is shown in Figure 1] and Figure 2. 
Note that the computer diagram is arranged such that the generator time response to changes indesign 


parameters can be quickly determined. 
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Figure 1: Response of Separately Excited Generator to Application of Exciter Voltage. 
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Figure 2: Response of Separately Excited Generator to Changes in Exciter Voltage. 
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TIME RESPONSE OF A TWO-WINDING TRANSFORMER 


Problem Description: 


To determine the time response of a particular design of a two- winding transformer for changesin 
the resistive load and input voltage. 





Equations : 
=e *2 
Miss 
C Ri 
IDEAL 
TR FOR 
SCHEMATIC ntiatiaiaiaal 
OF TWO-WINDING 
TRANSFORMER SIMPLIFIED EQUIVALENT CIRCUIT OF TWO-WINDING 
TRANSFORMER 
From the equivalent circuit the following primary and secondary voltage -equilibrium equations can 
be written : 
9 2 dI, 
By - (Ry + @ Ro) Ty ss (Ly + 2 ig) a Ey 
t 
Som om = ToR, 
By Ny 
— = — =a 
KE 6ONg 
Noly = Niky 
where 
Ry = Resistance of primary winding 
Ro = Resistance of secondary winding 
Ly = Inductance of primary winding 
Lo = Inductance of secondary winding 
a = Transformer turns ratio 
E;,k ..E, = RMS voltages 
Ij, I, = RMS currents 
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Computer Diagram: 


4 lO( Ly+a2Lo) Es 
+10 ! 











Results of Computer Solution: 
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Figure 1: Response of Output Voltage of Two-winding Transformer for Changes in Input Voltage. 



































Figure 2: Response of Output Voltage of Two-winding Transformer for Changes in Resistive Load, 
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INVESTIGATION OF UNSTEADY-STATE HEAT CONDUCTION 


Discussion: 


Many engineering problems involve the sudden heating or quenching of large slabs of material. In 


the annealing or heat treating of metals, plastics, glass, and rubber, for example, an accurate knowledge 


of the temperature vs. time and thickness relationships is required. Consider the following representative 


problem. 


Problem Description: ‘ 


Fal 





Determine the temperature distribution along the X 
direction of a long wide slab. The slab is initially at a 
high temperature (T]1) and is suddenly quenched in a 
large ice-water bath at 0°C. Assume that the temperature 
in the slab varies in the X direction only. 





FACE Equations: 
FACE =——* With the assumption made above, the unsteady-state 
heat transfercan be expressed by the following equation: 
x OT k 0?T 
: i 2 
| ee ; The boundary conditions are, for 7 = T (x, t) 
x=0 T Ww, 0) = f° O<a<B 
| wae T (0,t) = 0 t >0 
Where T (Bb, t) = 0 2 = Ff 
T = temperature in degrees C. 
t = time 
k = thermal conductivity of material 
Cp = specific heat of material (ave) 
p = density of material (ave) 


Using a finite - difference approximation of the second derivative of temperature with respect to the 


space dimension x, the second central difference becomes ........ 








ie. + _ (7 Taf 
dx? = AX? may © gad! 
where the subscripts 7 + 1, 7, 7 — 1 are referred to the number of cross section of the slab. 
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| 1.0 
DIMENSION X ————> X=B 


The use of the finite-difference approach is due to the fact that the Analog Computer can con- 
tinuously integrate with respect to only a single independent variable at any one time. This approach is 


equivalent to dividing the slab into sections as shown by the sketch and solving for the temperature at 
each section. 


Using the finite -difference approach the Analog Computer must simultaneously solve the following 
ordinary differential equations 

















(L$ - ztkz) {Ey - (2) 70) | T_ fe} = 0 
oa ioee) {123 - tis] + 7a | 
[2-9 (Az) fen -() [4] + ony 
jt ats ee) fe -(4) 4 (ryt 


rl} dT) _ (5) ‘a (2) (241 4 oa T, (t) = 0 
10 dt Cpp Ax? 10 10 


Computer Diagram: 





Results of Computer Solution: 
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Figure 1: Temperature Change vs. Time At Several Points Along X Direction. 


Application Note #8 (continued) 
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INVESTIGATION OF THE PROCESS PARAMETERS AFFECTING THE CONTROL 
OF A STIRRED TANK REACTOR 


Discussion: 

The continuous stirred tank reactor is an essential part of many chemical processes. The problem 
below shows, in simplified form, how a model of a CSTR may be set up and studied. 
Problem Description: 


To investigate the effect of parameter changes on the time 
response of the output concentration of a Stirred Tank Reactor. 


Equations : 
The reaction taking place is 
A x B 
When the specific reaction rate is defined as 





Material Balance: 








Material in v C; 
Material out w Gs 
Material reacted k& VrGe 
., Material accumulated yp Ss 
Summing : 
aoe. a ee 
dt Ver Ver 
Heat Balance: 
Heat in vpCp T; 
Heat out vp Cp To 
Heat transferred hA*(T, — T*) 
Heat due to reaction —QOV rk Cy 
Heat accumulated aie Bp ase 
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Summing: 





eee T= A (T, — Te) —QkGo _ 2 Tt. 
dt Ver p Vr Cp p Cp Vr 

where 

T,, T;, T* = temperature 

Vv = flow rate of material 

Vr = reactor volume 

t = time 

h = heat transfer coefficient 

A* = effective area of heat transfer 

p = average density of material 

O = heat due to reaction 

Cp = average specific heat 

C; = input concentration 

Co = output concentration 


Scaled Equations: 


The variables are scaled as shown 
possible but do not exceed +10 volts. 


[10C; 
LOG 


T; 
=. | 
LT 
j £2] 
10 
T* 
10 


[10K] = 





below so that the computer problem voltages are as high as 


] [10K] 
] 


“9 


10 Kan + (1006) [ 


a L10K] [10C,] 
Lio] -(¢ -) (106 - G) b9e0 _ Font 
l paler hA r 7? O [10K] [10C,] v Ts 
(2-42 -D-Gone OE 
Mig x Ve Ty GEE 10 10 100p Cp [10] Vr J 10 


Computer Diagram: 


C; —[10C;] 








VrB 
__@_ 
T; 
100 te 
we O 
Vv 
Vrp 





, [10K] [10C,] 
[10] 


Note that the inputs to the integrators have been multiplied by a time scale factor 5: For B < 1 


the computer solution will be 8 times faster than the response of the real physical system. 


Results of Computer Solution: 
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Figure la: Response of Output Concentration for 
Disturbance to Input Concentration. (v/Vr = 0.5 ) 
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Figure 1b: Response of Output Concentration for 
Disturbance to Input Concentration (v/Vp = 0.167) 
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Figure Ic: Response of Output Temperature of Reactor for Disturbance to Temperature of Material 
Entering Reactor. (v/V7 = 0.167) 


Application Note #10 (continued) 


ASX FP FPL I CATT ION ST Ss 





NUMBER 11 


ANALOG COMPUTE FR 


INVESTIGATION OF THE OUTPUT RESPONSE OF A POSITIONAL SERVO SYSTEM 


Problem Description: 


To determine the output response of a positional Servo System to input disturbances for changes in 
the system parameters. 
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FUNCTIONAL BLOCK DIAGRAM OF A _ SERVOMECHANISM 
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Equations: 
Error Device 6 — 0, = € 
Equalizing Network Ey} «KE; Tp+l 
€ alp+1 


Servo Amplifier Eg = KoE 3 = Ko (E]— Eo) 
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Computer Diagram: 





Results of Computer Solution: 
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Figure la: Response of Servo for step input (k; = 0.47, T = 
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Figure 1b: Response of Servo for step input (ky = 0.87, T = 
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Figure Ic: Response of Servo for step input (k, = 0.87, T = 
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INVESTIGATION OF THE RESPONSE OF A NUCLEAR REACTOR 


Problem Description: 


To investigate the time response of a Design of a Nuclear Reactor for changes in the Reactor mul- 
tiplication factor dk. 


Equations : 


The equations! representing the performance of a Nuclear Reactor are as follows: 


_ i= 6 





. =F n+ % Ay C; 
dé ] * i= |] 
d Ce B 
a Ay C; 
dé }* 


To represent a Nuclear Reactor with five delay groups the following differential equations must be 
solved simultaneously : 


1=5 
dn = dk —B % + dj C; 
dé \* i= 
dC, 8B 
abel. = i: n— A, Cy 
, dé | * 
dC 
_ _—2 7 Be n — AgCo 
dé i# 
dC3  B3 
dé ] * 
dC, Ba 
— «= — nan rz Cy 
dé i* 
dCs Bs 
———S SS 6, Oe As Cs 
dé | * 





M.A. Schultz, ‘*Control of Nuclear Reactors and Power Plants’’, McGraw Hill Book Co., Inc., New York, 1955. 
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Results of Computer Solution 
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INVESTIGATION OF THE REACTION OF SODIUM VAPOR DIFFUSING INTO A HALIDE 


Discussion: 


Many chemical reactions are so rapid that diffusion of the reactants and products around the re- 
acting zone will completely control the system behavior. Such is the case when sodium vapor reacts 
with a halide. In the problem below, vaporized sodium was admitted through a nozzle into a large reactor 


filled with a halide gas. The equations describing this process demand a trial and error type method of 
solution. 


Problem Description: 


To investigate the effect of changes in design parameters on the reaction of Sodium Vapor 
diffusing into a Halide. 


Equations : 


The following equation! describes the reaction of Sodium Vapor diffusing into a Halide: 


dV  KV(V+x+x,—-1) 


d x? we + My 





The boundary conditions are 
V=1— x, at x = 0 
V = 0 as Xx > © 


The parameters (1 — x,) and K correspond to a measure of the nozzle opening used in a reactor 


and chemical rate of reaction, respectively. The solution to the equation is obtained when dV/dx and V 
become zero simultaneously. 


Scaled Equation: 


The independent variable x will be replaced by time so that x = ¢. The variables are scaled as 





follows: 
og vl 0 < & < + 1 
dt? 
190 2h 1. ¢#, < + W 
dt 
[10V ] 


The scaled voltages do not exceed the maximum permissable computer voltages, +10 volts, and 
permit the following scaled equation to be written: - 


1. ‘*Mathematic Tables and Other Aids ta Computation’’, Physics Abstracts, Vol. IX. July, 1955, pp. 112-116 
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Computer Diagram: 
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Vo 7 (1-X)) & -|0 | +X 
lo(i~x,) . 
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sie 
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For 


Note that the inputs to the integrators have been multiplied by a time scale factor 
8 > 1 the computer solution will be B times siower than the real reaction of the sodium vapor. 


Results of Computer Solution: 
The results of several trial runs are shown in Figure 1. From these trial solutions it can be 
seen that if V, is too large, V will change sign before V. If Vo is too small V changes sign before V. 


The objective in choosing a Vp» is to cause both V and V to go to zero simultaneously. The best tech- 
nique to use is to find a “‘small’’ value of V,> and then a “‘large’”’ value of Vo, and approach the “‘small’’ 


value by successively reducing Vp until the correct solution is obtained. 


Typical final solutions for V as a function of the independent variable are shown in Figure 2. 
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Figure 1: Typical run of trial solution for K = 0.1, 1 — x, = 0.20. 
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Figure 2: Typical Results of Calculations of the Reaction of Sodium Vapor diffusing into a Halide. 
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